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1. Introduction 

With the widespread availabihty of powerful lasers, the interaction of atomic systems 
with intense electromagnetic radiation has become the subject of much study. A 
variety of multiphoton phenomena have been observed, remarkable by their non- 
linearity. The field may induce new processes that can only occur through the 
absorption or emission of more than one photon. Multiphoton ionization, above- 
threshold ionization and harmonic generation have attracted most interest, while 
electron-atom collisions in the presence of a laser field have received comparatively 
less attention. These are commonly described as laser-assisted, in that the field is not 
required for scattering to occur, although its presence can give rise to new inelastic 
processes such as simultaneous electron-photon excitation. 

Relatively little experimental work has been devoted to electron-atom collisions 
in a laser field. A number of groups have investigated clastic or free-free scattering 
on rare gas targets (usually argon) in a CO2 laser field at intensities up to 10^ 
W/cm2 (Andrick and Langhans 1976, 1978; Weingartshofer et al 1977, 1979, 1983; 
Bader 1986; Wallbank et al 1987a, b, 1992). Measurements have been performed 
over a wide range of incident electron energies, mostly at large scattering angles 
corresponding to backward scattering, for geometries in which the polarization axis is 
either parallel to the momentum of the incident electron, or parallel to the direction of 
momentum transfer. More recent experiments (Wallbank and Holmes 1993, 1994a, b) 
measured angular distributions for geometries in which the polarization axis is almost 
perpendicular to the momentum transfer. The first experimental investigation of 
simultaneous electron-photon excitation was performed by Mason and Newell (1987, 
1989) on helium in the field of a CW CO2 laser at intensities from 10^ to 10^ W/cm^ 
This was soon followed by studies at intensities up to 10^ W/cm^ (Wallbank et al 
1988, 1990). The only experiment at a higher frequency is that of Luan et al (1991) 
who used a Nd-YAG laser of intensity 10^*^ W/cm^. Details can be found in the review 
by Mason (1993). 

A larger number of theoretical studies have been undertaken, and are 
comprehensively reviewed by Ehlotzky et al (1998). The challenge for theory lies in 
accurately treating both radiative couplings and the dynamics of the collisional system 
over very large distances: each of the electron-target, laser-electron and laser-target 
interactions needs to be taken into account. If one of these dominates, a perturbative 
approach can be adopted. In early work on fr(X'-fr(x; scattering for example, Bunkin 
and Fedorov (1966) treated the laser-electron interaction to all orders using Gordan- 
Volkov wave functions, and the electron-target interaction to first order using the 
Born approximation. The dressing of the target by the field was neglected. KroU 
and Watson (1973) included higher order terms in the Born series to obtain a formula 
valid when the scattering potential is weak or when the frequency of the laser field 
is small (the low frequency approximation). An alternative approach for treating 
the laser-atom interaction non-perturbatively is furnished by Floquet theory. This 
is applicable to periodic Hamiltonians and reduces the problem to be solved to that 
of a set of time-independent coupled equations. The Floquet expansion has been 
combined for example with a partial wave expansion for potential scattering (Dimou 
and Faisal 1987), and with the Lippmann-Schwinger integral equation to study low- 
energy free-free scattering by various potentials (Kylstra and Joachain 1998, 1999). 
Most work on inelastic processes has concentrated on high collision energies using 
methods based on the Born approximation (Rahman and Faisal 1976, Jetzke et al 
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1984, 1987, Francken et al 1988). Lower energies have been considered mainly within 
the low-frequency approximation, which is valid for CO2 lasers (Gcltman and Maquct 
1989, Mittleman 1993, Fainstein et al 1995) but not for Nd-YAG. An attempt to 
go beyond this approximation using a pseudo-potential model has been discussed by 
Faisal (1987). 

There now exists a variety of powerful computational methods for the study 
of electron collisions with complex atoms. i?-matrix theory for example has been 
developed over many years (Burke and Bcrrington 1993, Burke et al 1994) and is 
now used regularly for the production of vital atomic data needed in a variety of 
applications. The main feature of this theory is the division of configuration space 
into inner and outer regions, in each of which the wave function can be expanded 
in a locally adapted basis. The global wave function is then built by matching the 
solutions at their common boundary, using the inverse of the log-derivative matrix 
or i?-matrix. An extension of the method to multiphoton processes was proposed by 
Burke et al (1990, 1991, hereafter referred to as paper I), using the Floquet expansion 
to describe the interaction of the atomic system with a laser field. A practical 
solution in the outer region was developed by Dorr et al (1992), hereafter referred 
to as paper II. Subsequently, this i?-matrix Floquet (RMF) theory has been used to 
study multiphoton ionization of numerous atoms and ions, laser-induced degenerate 
states (LIDS) and harmonic generation (see for example the review by Joachain et al 
2000). A version describing multiphoton processes involving diatomic molecules has 
recently been developed (Burke et al 2000) and applied to multiphoton ionization of 
H2 (Colgan et al 2001). 

In papers I and II, the RMF theory was presented mainly with multiphoton 
ionization in mind. Its application to laser-assisted scattering has been described in 
detail only for the case of a potential (electron-proton scattering, Dorr et al 1995). We 
have since extended the method to laser-assisted scattering by atomic targets. We first 
applied it to the study of free-free processes in e~-H (Charlo et al 1998) and e~-He 
(Charlo 1999) collisions in a CO2 laser field. We showed that the differential cross 
sections display a deep minimum when the momentum transfer is nearly perpendicular 
to the direction of the laser field, in accordance with the low frequency approximation. 
More recently, we have studied the simultaneous electron-photon excitation of helium 
in a Nd-YAG laser field (Terao-Dunseath et al 2001), for which the low frequency 
approximation is no longer valid. We showed that the cross section is dominated by 
the He~(ls2s^ ^S) resonance and that a strong AC Stark mixing occurs between the 
ls2s "^S and ls2p "^P" states. The purpose of this paper is to present the modifications 
that needed to be introduced in the original i?-matrix Floquet theory in order to 
identify unambigously field-dressed target states and to calculate cross sections for 
laser-assisted scattering by a general atomic or ionic system. 

Atomic units are used throughout. We shall also use the Fano-Racah phase 
convention for the spherical harmonics, with 

where CS implies the usual Condon-Shortley phase. 

2. The Floquet approach: basic equations and properties 

We consider the scattering of an electron by an atomic system composed of an infinitely 
heavy nucleus of charge Z and N electrons in a linearly polarized laser field. We 



RMF theory of laser-assisted electron-atom scattering 



4 



suppose that the duration of the laser pulse is much longer than the collision time 
and that the field is homogeneous over the interaction region. We therefore adopt the 
dipole approximation with a field described by the vector potential 

A{t) = z Ao cosujt (1) 

where u is the angular frequency and we have chosen the Z-axis parallel to the direction 
of polarization. 

The wave function for this system satisfies the Schrodinger equation 



i^*(X^+i,i) 



i N+1 , ' 

- -A{t) • ^ Ve + -i^A\t) 

e— 1 



■^{XN+i,t) (2) 



where X n+i = {xi, X2, ■ ■ ■ ,xn, a^Ar+i} denotes the set of space- and spin-coordinates 
of the -I- 1 electrons with = {re,cre}, and Hm+i is the field- free Hamiltonian 
operator 



AT -1-1 

Hn+1 = ^ 

e=l 



-V 

2 " 



N+l 



y 



e'>e=l 



The presence of the linearly polarized laser field breaks the spherical symmetry 
of the system by introducing a preferred direction in space so that the total atomic 
angular momentum L is no longer defined. Its projection M along the polarization 
axis is however conserved since the system is invariant under rotation around this 
axis. The dipole operator connects atomic states of opposite parities, so that equation 
(||) separates for well-defined values of the total atom-laser parity 11. As we neglect 
relativistic corrections, the total atomic spin S and projection Mg are also good 
quantum numbers. In what follows, we shall work in subspaces with fixed quantum 
numbers A/nS'Afg, which will usually not be explicitly indicated. 

As the Schrodinger equation (^) includes a time-dependent potential which is 
periodic, it admits solutions in the form of a Floquet-Fourier expansion 

oo 

*(Xjv+i,0=e-'^* ^ e-'"'^*vi,„(XA.+i) (3) 

n— — oo 

where E is the quasi-energy of the solution. The conservation of total parity 11 implies 
that each Floquet component in the expansion (H) includes atomic wave functions of 
well-defined parity, alternating with the parity of the Floquet index n. By substituting 
(||) into (H) and projecting on a particular Floquet component, we obtain an infinite 
set of time- independent coupled equations for the ^'„(Xjv+i) that can be solved using 
standard methods. 

Another general property of the Floquet-Fourier type solutions (||) is the Shirley 
symmetry (Shirley 1965): 

oo oo 
_ Q-iEt 'y ^ ^-inujt^E _ ^-i{E+mLij)t 'y ^ ^-i{n-m)ujt^^E 

n— — oo n— — oo 

oo oo 
^—i{E+muj)t ^ ^ ^—inuty^E ^—i{E+niuj)t ^ ^ ^ — inuty^E+muj 



e 

n— — oc n— — oo 



This implies that equivalent sets of solutions of the Schrodinger equation (||) exist, 
each characterised by quasi-energies with period uj and Floquet components identical 
apart from a shift in their indices. 
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3. i?-matrix inner region 

In i?-niatrix tlieory, the inner region is defined as the portion of configuration space 
where the radial coordinate of each electron is smaller than the radius a of the sphere 
encompassing the target states. The treatment of the problem in this region, which 
is analogous to that of atomic bound states, has been fully presented in paper I. 
It is however useful to recall the equations that are essential for understanding the 
treatment of laser-assisted scattering in the outer region. 

Since all the radial coordinates are small, it is appropriate to define the wave 
fimction in the length gauge by the unitary transformation 

/ . N+l \ 

^{XN+i,t) ^ exp i-'-A{t) ■ ^^{XN+i,t). (5) 

\E''"(XAr_|_i, t) is solution of the Schrodinger equation 
d 



i^*L(XAr+i,<) 



N+l 

HN+i+£{t)- Vre 



■^''{XN+i,t) = H^-^^{XN+i,t) (6) 



e=l 

where the electric field £ is defined as 

c = -— = zt{) smujt (J) 

c at 

with £q — wAo/c. By substituting the Floquet-Fourier expansion of ^'^(Xjv+i, t) 

oo 

vI/L(XA,+i,t) = e-'^'* ^ e-"'^'^^^{XN+i) (8) 

n— — oo 

into equation (^, multiplying by a particular exp(iria;t) and integrating in time over 
a period 27r/cj, we find that the Floquet components '^\-^{Xm+i) are solutions of the 
coupled differential equations 

(Fjv+i -E^- nuj)^^(XN+i) + Dn+i^\,+i{Xn+i) - Dn+i^\,^i{X n+i) = 
with 

£o 



N+l 



Dn+1 ■ '^e- 

e— 1 

We first determine a discrete i?-matrix basis {^'^(Xjv+i, t)} satisfying the 
usual i?-matrix boundary conditions. This is done by diagonalizing the hermitian 
Hamiltonian TJp -|-Lb, where Lb is a Bloch operator cancelling the surface terms at r = 
a arising from the kinetic operators. The Floquet-Fourier components '^\^^^{X n+i) 
of the i?-matrix basis functions are expanded in a set of fully antisymmetrised 
(A'^-f l)-electron functions which are formed by coupling channel functions with defined 
symmetry F — LMttSMs, (where tt is the parity of the (N-l-l)-electron system) 

^uIVn+i)^ Y,{UMdm\LM){S^Ms,\ii\SMs) 

Mim Ms lJ- 

X <^i{XN)Yi^{eN+l,(pN+l)x^f,{o-N+l) (9) 

to continuum orbitals UneifN+i) which are chosen as regular eigenfunctions of a field- 
free Hamiltonian including a static potential of the target ground state, with zero log- 
derivative at rjv+i = a. These orbitals are furthermore built to be orthogonal to those 
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in the target states ^i{X]y). {N -\- l)-electron bound configurations constructed using 
the target orbitals are also included in order to account for short-range correlation. 

By projecting (^) onto the channel functions (^) and onto the n"* component in 
photon space, and evaluating at the boundary tat+i = a, we obtain the relation 

v' 

(where v represents the set of indices {Ti£n}) between the amplitudes of the radial 
functions 

i^, = ($i;|5(r^+i-ra)vI/L) 
and their derivatives 

i^, = ($i;|5(r^+i-ra)-^vI,L) 

OTN+l 

in terms of the i?-matrix elements 

lyT^^fcT^ (10) 

which depend on the quasi-energies E^ and the surface amplitudes of the eigenstates 

7.fe = ($j;i<S(rAr+i-ra)vI/L). (11) 

In principle the summation in ( p^ is infinite, in practice it is truncated. An 
approximation for the truncated part is provided by a Buttle correction (Buttle 1967). 



4. Close-coupling equations using field-dressed target states 

The outer region is defined as the portion of configuration space where the N electrons 
of the target have a radial coordinate smaller than a while that of the coUisional 
electron is larger. Exchange between the coUisional electron and the target electrons 
is negligible. From now on, we can drop the index A'^ + 1 from the coordinates of the 
coUisional electron without possibility of confusion. 

To solve the Schrodinger equation in the outer region, we adopt a close-coupling 
expansion of the total wave function. Using the same channel functions (|^) as in the 
inner region, we could write 

oo 

vI/L(XA.+i,<) = e-'^'* e-"-* 5] $r (V)-FrH,„(r) (12) 



r 



ri= — oo Livii 

where the notation \r denotes the coordinates X^v of all the target electrons together 
with the angular and spin coordinates of the coUisional electron. Substituting ( p^ 
into (^, we find that the radial functions F^-^^{r) satisfy the set of coupled ordinary 
differential equations in matrix notation 

i;-^ + ^ + .^>^W(.)F (13) 
dr^ r"^ r / 

where z = Z — N is the residual charge of the target. All quantities operating 
on F on the left-hand side of ( p^ ) are diagonal matrices. The elements of are 
the square of the channel momenta. The matrix W(r) includes the usual long- 
range multipole potentials a^^-'/r^ arising from the interaction of the coUisional 
electron with the target, together with the couplings of all the electrons with the 
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laser field. The interaction of the target with the field does not depend on the 
coordinate of the coUisional electron and therefore gives rise to non-diagonal constant 
couplings in equation (p^. In paper I (appendix A2.1.2), these terms were removed 
by diagonalizing the constant coupling matrix, which amounts to defining channel 
functions with the field-dressed target states coupled to the angular and spin wave 
functions of the coUisional electron. This approach has two disavantages. Couplings 
in the outer region are first defined in the channel basis (|9|) where their expression 
is cumbersome, then transformed into the field-dressed target basis. Secondly, target 
states differing only by the sign of their magnetic quantum number Mi are degenerate 
so that the diagonalization procedure can lead to any linear combination of the 
solutions with ±|Mi|. 

We now show that the equations in the outer region are greatly simplified by first 
determining the field-dressed target states and using these to construct the channel 
functions. The magnetic quantum numbers of the target and therefore of the coUisional 
electron are defined unambigously and scattering cross sections can be obtained in a 
straightforward way. 

The field-dressed target wave functions \l'^(Xjv,i) satisfy the time-dependent 
Schrodinger equation (0), with -I- 1 replaced by N: 



i^*^(Xw,i) = 



N 

HN+e{t)-"^re 

e=l 



^^iXN,t). (14) 



We adopt the Floquet-Fourier expansion 

oo 

^^{XN,t)=e-'^^* e-'""*<(XAr) (15) 

n— — OG 

and expand each '^'^{Xn) over the eigenfunctions ^i{XN) of Hn used in (^: 

vl/^(X^)-^<i>,(X^)af„. (16) 

i 

The expansion (Hq) includes atomic states with different Li but the same Mi, in, 



Si and Msi- Substituting (|15|) and (|16| ) in (|14D and projecting on a particular 
ex-p{—inuit)^i{X n) , we find that the coefficients af„ are solutions of the eigenvalue 
problem 

(w, ^E^ - nuj)al, + ^($, | Dn \ $»')(aL+i " aL-i) = 0, 

i' 

where we have used HN^i{X n) — Wi^i{X f^). As a consequence of the Shirley 
symmetry (^, the quasi-energies can be grouped in sequences of period lu. The 
presence of the linearly polarized field breaks the degeneracy in the magnetic quantum 
number Mr but states with the same value of \Mt\ are still degenerate. Since the 
Hamiltonian (|lj) is non-relativistic, states differing only by Mg^. are degenerate. 

We define channel functions with well-defined quantum numbers M, 11, S and Ms 
by coupling the 'i''^{XN,t) to angular and spin functions of the coUisional electron: 

§T^™(\r,i) = e-'^-* J2 e-""* ^ (^tM^, I 5A/5X(X^)>£„(^?, ^)xi ^(a). 

(17) 



-C50 
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In this basis, the close-couphng expansion of the wave function ^^{X ]y^i,t) is 

vI/L(XA,+i,t) = e-'^'*5]e'^-*5]$T.™(\r,t)iFL,,„(r). (18) 

T tm 

By comparing equations ( p^ and (|lj), we find that the radial functions F^^.^{r) are 
related to those in ( p^ ) by the transformation 

Frwuir) = ^Ortt'n,ThnFTim{'r) 

Tl. 

with 

Orit.'n,Ttm = {LiMilm\LM)aJn5i'i5M,MT- 

The Kronecker symbol (5^'£ arises since this basis transformation concerns only the 
target and not the coUisional electron, while the second reflects the fact that the 
magnetic quantum number of the target is not affected. If we denote by F the vector 
of i^pj£„ and by F the vector of F}ji^^^ we have in matrix notation 

F = OF 

and similarly for their derivatives 
F = OF. 

Since the i?-matrix ( [lO| ) corresponding to the channels (||) relates F and F by 
F = RF, 

the i?-matrix R corresponding to the field-dressed channels ( p7| ) is related to R by 
the transformation 

R = O'RO 

where O* is the transpose of O. The long-range multipole potential coefficients in 
equation ([l^) must also be transformed in a similar way: 

5(A) ^0*a(^)0. 

The channel energies k? /2 must be recalculated in this new basis, while the centrifugal 
term remains the same. 



5. Length to velocity gauge transformation and propagation in the 
velocity gauge 

In the outer region, the potential in the length gauge S{t) ■ r diverges as the radial 
coordinate of the coUisional electron increases. It is therefore more appropriate to 
express the interaction of this electron with the laser field in the velocity gauge, where 
the potential remains finite. We conveniently perform this gauge transformation at 
r = a, although in principle it could be done at a larger distance. The relation 
between the wave function in the length gauge ^'^(Xjv+i, t) and that in the velocity 
gauge ^'^(XAr+i,t) is 

*^(Xw+i,t) =exp(-^ f A\t')dt' ~-A{t)-r\^'^{XN+i:t). (19) 
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The Schrodinger equation in the velocity gauge is 

JV 

Hn + £{t) ■ 



i|-*V(X^.+i,i) 



dt 



e=l 
N 



e=l 



r ^ — ^ I 7* — 7* I r 



^^{XN+i,t). (20) 



Analogous to (18), we expand "if^ {X ]y+i,t) as 

vI/^(X^.+i,i) =c-'^"*5]e'^-*5]$T.™(\r,i)-^^^,™(r-). (21) 



hn 



Substituting ( |2l| ) into equation ( |20| ) and projecting on a particular channel function 
(p^, we obtain the close-coupling equations in the velocity gauge. In matrix form, 
they are 



d' £{£ 
dr^ r 



2z 



K 



Ar 



A„ex ~(A) 



(22) 



A=l 



where Pd/dr and or ' jr arise from the term Aif) ■ V. The elements of P are 

^0 
C 



'(£'-TO)(^' + m) 
(2£' + l)(2f -1)' 



' (f - m)(£ + m) 
(2£+ 1)(2£- 1)' 



while those of a'-"'^-' are 



^0 



mm 



' {£' - m){e' + m) 
(2f + 1)(2£'-1)' 



' (^- m)(^ + m) 
(2£+ 1)(2£- 1)' 



The Kronecker symbol St,t'±uj means that T and T' are the same field-dressed target 
states (i.e. with the same atomic structure), differing only by their quasi-energies 
Et — Et' ± oj- The dipole operator couples channels where the collisional electron 
absorbs or emits one photon. P and a^^^ obey the usual one-electron dipole selection 
rules A£ = ±1 and Am = 0. The fact that P and 5*-^^ are real is a consequence of 
the choice of phase for the field and the Fano-Racah phase convention adopted for the 
spherical harmonics. Since the Hamiltonian is hermitian, P is antisymmetric and a^^-* 
is symmetric. 

In order to solve the close-coupling equations (|2^) we must first transform R into 
the velocity gauge. Using equations ([l]) and (|l9|), we write 



*^(Xjv+i,i) = exp<ji^tlexp<|i-^sin2wt- -AqZ • rcoswtl *^(^Ar+i,i) 



4c2 



a: 



(23) 



By comparing this with equation (^ ) and using (||), we may first identify 



E'' = E^-^ = E'^-Ep 
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where Ep is the ponderomotive energy. The second exponential in equation 
periodic in time and can be written as the product of two Fourier series, 



exp < 1 



A, 



Suic' 



° sm2ujt AQZ-rcosiot 



/ 42 
/ ^0 



E 



^0 



(24) 



Substituting equatio n (p^ ) into (2^, replacing and by their respective close- 
coupling expansions (|18[) and (pT) and projecting onto the channel functions (|l^), we 
obtain relations between the radial functions F^gmi^) ^remi^)- These can be 
written in matrix form as 



= AF 

where the elements of A are given by 

/ 42 



(25) 



■t'-t-k-2J\f 7 



AT 



(2£ + l)(2f + 1){£ - m)!(f - m')l 



7r/2 



{£ + m)\{e' + m'y. 
Ao 



e] <5„ 



d6isin6iP7'(cos6l)P;r (cos6') J_fc_2A^ 

\ c J 

(26) 

Since the gauge transformation applies only to the scattered electron, the channels 
coupled by the matrix A are those where the atomic structure of the field-dressed 
target state is the same: Et = Et' + few, where k is an integer. The Kronecker 
symbol Emm' arises from the fact that the laser field is polarized along the Z^-axis. 
The integral in equation (^^ is non-zero only ii £-\-£' — k — 2J\f is even, implying that 
the matrix A is real. 

Since the radial functions and F are related by the simple linear 

transformation (p5[), it is easy to show that 

^ A fA + AR-i 



where the elements of A are the derivatives of ( |26| ) at r = a. In order to apply standard 
propagation techniques to solve (|2^), we remove the Pd/dr terms by transforming the 
radial functions as 



-Pr/2Q^ 



(27) 



This is analogous to a Smith-type diabatization (Smith 1969) in molecular scattering. 
Substituting ( p7[ ) into (p2[), we find that the functions G satisfy the set of coupled 
equations 



d^ 
dr2 



,Pr/2 



2z 



A^. ~(A) 

E 

A=l 



-Pr/2 







where the couplings oscillate with the distance r. The i?-matrix corresponding to the 
functions G is given by 

-1 



,-Pr/2 
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which is symmetric. A detailed discussion of the properties of has been presented 
in paper II. 

To propagate the matrix R^j, we first inverse it to obtain the log-derivative 
matrix. We propagate this outwards using the method of Johnson and Manolopoulos 
(Manolopoulos 1986) which allows oscillating potentials to be treated with relatively 
large propagation steps. The propagated log-derivative matrix is then inverted back 
to the matrix Rg. 



6. Scattering boundary conditions and asymptotic solutions 

After propagating the i?-matrix in the velocity gauge, we need to match to the 
asymptotic solutions corresponding to the scattering boundary conditions. These can 
easily be defined if the close-coupling equations are asymptotically uncoupled, which is 
not the case in the velocity gauge. We apply the Kramers-Henneberger transformation 
(Kramers 1956, Henneberger 1968) 

^^iXN+i,t) - exp[-ia(i) • p] ^f^iXN+ut) = ^^{Xn, r + a{t), tr, t) (28) 
where 

1 /■* 

alt) — - I A[t')(lt' — z aosmwt — z Eo/lo^ smuot 
c J 

is the oscillation vector of a free electron in the laser field. It must be emphasized that 
this transformation is applied only to the coUisional electron. The target electrons, 
which are tightly bound to the nucleus, are kept in the length gauge. By substituting 
~8|) into (po|), we obtain the close-coupling equations in the acceleration frame: 

N 



e=l 

i_, z - 1 



^^{XN+i,t). (29) 



2 \r + a{t)\ ^|r + a(i)-r 

In the limit r —^ oo, ag and are negligible with respect to r. The right-hand side of 
( p9| ) reduces to the target Floquet Hamiltonian plus the Hamiltonian of an electron 
in the screened Coulomb potential of residual charge z — Z — N. Two sets of linearly 
independent solutions can therefore be defined, with boundary conditions 

^tiXN+ut) ^ e-'^'* V y2iSTMs^^^i\SMs)e'''-'^^iXN,t) 

r—>oo ' ^ ' ^ 2 

1 ie±'^-Wr,™((?,^)xiJ'^) 



which can be rewritten as 

^>t{XN+i,t) — > V V(5tMs i/i|5Ms)f^(X^,t) 



xe -.^*^ie±'''-'Wy,„(0,^)xi^(^) (30) 



where e,^ — — Et — k't/'^, is the channel energy of the scattered electron in the 



acceleration frame. From equation ( ^81 ) and the Floquet expansions of \E'^(Xjv+i, t) 
and \l/^(Xjv+i, <), it is obvious that the quasi-energy in the acceleration frame E^ 
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is equal to . We have introduced the coUective index ly = {T,£,m} to identify 
the channel. The Coulomb phase is Oi,{r) — k^r — ^7r/2 — 7/^ h\{2k^r) + crpX^v), with 
= -z/k„ and (t<i(?7^) = argr(£ + 1 + ir]^). 
We now define the boundary conditions in the velocity gauge by applying the 
frame transformation (ESh to (BOf): 



^^{XN+i,t) Y.T.^^TMs^^^^\SMs)^^{XN,t) 



1 1 



X e---*^-e±'«"('-°)y,™(0„,^) xy{<T) (31) 

where (t-q, 6*0,, (/j) are the spherical coordinates of the shifted vector ra{t) — r + a.{t). 
In the limit r 00, we have 



7'q = r — ao sin cut cos + O 

so that 



1 = 1+0^^ 



cos 6*0, = cos -\- O 



The azimuthal angle of is simply given by 
We can therefore use the approximation 

expjifcrQ — iriln{2kra)} — expjifcr — i77ln(2fcr)} exp{— ifcao sinwt cos0} + O (^~^ ■ 

The time-dependence in the phase 0^{ra)leads to a mixing of the channels in the 
velocity gauge. We may write equation {pW as 

^^AXN+i,t) ^ ^$,,(V,t)A(?)e--^*^ie±'«"W (32) 



where v' = {T',i' ,m'}. By equating (|3l|) and (32) and using 



J„(x) ^ — / d6'e'(^'''"^-"*) 
27r Jo 

we find that 

a'^°1 = fdn Yf1^,{0,ip) Y,^{0,ip) Jfe(Tfc.aocos0) 5,, 



where k, defined by Et' = Et + kuj, must be an integer if A^^}^ is to be non-zero. 
This is again a consequence of the fact that the frame transformation is applied only 
to the scattered electron. The magnetic quantum numbers are conserved since the 
transformation is along the Z-direction. 

We now use (|3^ ) as the first term of an asymptotic expansion of the radial 
functions in the velocity gauge 

*y (X^.+l, t) = Yl ^'^'(V^ E -A'l e-'-'^- e±'^^ W. (33) 



" " ^fk r 
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In matrix notation, we have for the radial functions 

fV(^)^ 1 IgiieM y 1a('^) (34) 
where the factor in front of the summation denotes a diagonal matrix with elements 



exp(±i6'jy(r))/vAJ^r. By substituting ( p4| ) into (22), we obtain a set of recurrence 



relations for A*^''^. Details have been given in paper II. 
7. Scattering amplitudes and cross sections 

The close-coupling equations in the acceleration frame are asymptotically uncoupled, 
as in the field-free case. We can therefore impose the usual scattering boundary 
conditions 

F^i^'ir) — > fci7 (sin 6*^(5^1,/ -I- cos 6'i. if 1.^/), fc^ > 

r—)-oo 
r — ^oo 

The if-matrix is calculated by matching the solutions propagated from the inner 
region to the asymptotic solutions in the velocity gauge. The matching equation is 
obtained by writing the propagated solutions as 

n-t-Tio 
v" = l 

where the elements of H are the asymptotic solutions with boundary conditions in the 
acceleration gauge 

Huy'{r) > sinO^ 5^1,1 i/ = 1, . . . n; i/' = 1, . . . no 

r — >oo 

Huu'{r) — > C0S6',, i/ = 1, . . . n; i^' = rio -I- 1, . . . 2no 

Hvu'{r) > Suu'-n^ 1/ = 1, . . . n; I/' = 2no 1, . . .71 -I- no. 

r — 'oo 

The first two sets correspond respectively to the imaginary and real parts of the 



solutions (33) with kf, > 0. For the open channels, we need to calculate only one 
of the exp(±ifc^r) solutions. The third set corresponds to the solutions ( ^3|) with 
< 0. For the closed channels, k^ is imaginary and we calculate the solutions 
decaying asymptotically. The matrix H is of dimension n x (n + no), where n is the 
total number of channels and no is the number of open channels. The if-matrix is 
obtained from the solution X to the equation 

MiX = M2 

where M2 and Mi correspond respectively to the first no and the last n columns of 
the matrix 

P 

H — Rg H — Rg— H. 

The matrix X has dimension n x Hq and the if-matrix is given by the first rows of 
this matrix. 

The scattering amplitudes are obtained from the if-matrix in a manner similar to 
the field-free case (see for example Smith 1971). Differences arise due to the presence 
of the laser field which breaks the spherical symmetry. We first suppose that the target 
is neutral and that the incoming electron has a momentum ki and a spin projection 
/Xi while the outgoing electron has a momentum kf and a spin projection fif. The 
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scattering amplitude for a transition from the initial field-dressed target state i to the 
final field-dressed target state / is found to be 

^Jp^Y^iS^Ms^l^, I SMs){SfMs,l^^f I SMs) 

xT.J2 Y;^{k)Yi>m'{ks) T.'fl^^,. (35) 
mh w 

'^iimfi'm' element of the T- matrix which is related to the 5- matrix by T = 1 — S, 
with S = (l-|-iK)(l — iK)^^. The T-matrix is block-diagonal with respect to the three 
good quantum numbers Af , 11 and S. It is independent of Ms and the sign of M. 
The expression (^5|) is simpler than that in the field-free case since the total angular 
momentum L is no longer a good quantum number. It is therefore not necessary 
to uncouple the angular momenta of the target and scattered electron to define the 
scattering amplitude. 

The differential cross section for the scattering of an unpolarized electron by an 
unpolarized target is obtained by summing the square of the scattering amplitudes 
over all possible final states and averaging over all initial states: 



^ .2 



25, + 1 ^ 



'^(25+1) 



(36) 

A/n w 

where Ei is the kinetic energy of the incoming electron. The total cross section is 
obtained by integrating (|3^) over the outgoing angles 

* * S MUU'l" 

where we have m — M — Mi and m' — M — Mf. 

In the case of an ionic target with residual charge z, the scattering amplitude 
contains two terms 

hti.,ftif{K,kf) = fc{k„kf)5^f + f['^^j^^{k,,kf). (38) 

The first term is the Coulomb scattering amplitude 

fc(k,,kf) = ^ g2i.o(/c)gi(./fc)l„[sin^(e/2)] 

■'^^ " 2fc2sin2(e/2) 

where Q is the angle between ki and kf and k is the amplitude of ki and kf. The 
second term is similar to (K5f) and has the form 



ft:j,,ik,,kf) ^-^Y.{S,Ms^^fM I SMs){SfMs,l^^f I SMs) 

X E E e'[^'('=-'+^"('^^)ir/„(fe,)r,,„,(fe/)7;tL"f.'™'- (39) 

The differential cross section is obtained by taking the modulus square of (p8|), 
averaging over the initial states and summing over the final states. As in the field-free 
case, the total elastic cross section is not defined as it diverges. 

It should be stressed that the above formulae are for transitions between field- 
dressed target states. The Shirley symmetry implies that these can be grouped in 
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equivalence classes, each with the same atomic structure and a quasi-energy spectrum 
of period oj. When the AC Stark mixing is weak, these quasi-energies will lie close to 
the field-free thresholds shifted by an integer multiple of the photon energy. It is then 
natural, though not necessary, to use these integers as labels, with the field-dressed 
states closest to the field-free ones defining the "zero-photon" thresholds. A transition 
from state (a, M) to state (6, N) can therefore be interpreted as the atomic transition 
a b with net exchange of \AI — N\ photons. The coUisional electron emerges with 
a kinetic energy differing by \M — N\uj from the kinetic energy associated with the 
zero-photon transition. When the AC Stark mixing is very strong however, it may 
be difficult to identify the field-dressed threshold as a field-free threshold shifted by 
a particular number of photon energies. For inelastic processes, we may choose a 
particular transition between two field-dressed states i and / as the "zero-photon" 
transition for the associated atomic states. This can then be used as a reference for 
counting the number of photons exchanged for transitions between other members of 
the two sequences containing i and /. 

8. Computational implementation 

An advantage of the i?-matrix Floquet theory is that it can re-use existing _R-matrix 
inner region computer codes for field-free electron-atom collisions and photoionization 
that have been developed over many years (Berrington et al 1978, 1987). The Fano- 
Racah phase convention for the spherical harmonics used in those packages has to be 
kept, which yields rather unusual expressions for the dipole matrix elements. This 
also justifies the choice of phase for the laser field (^, since a real symmetric Floquet 
Hamiltonian is more convenient from the computational point of view. 

The first task is to define a set of field-free target wave functions. While in 
principle any method could be used to determine these, in practice only the structure 
programs CIVS (Hibbert 1975) and Superstructure (Eissner et al 1974) are interfaced 
with the i?-matrix codes. Once the target state energies, wave functions and dipole 
matrices have been obtained, the field-free i?-matrix inner region codes are used to 
calculate dipoles, energies and surface amplitudes for the (N-l-l)-electron system. The 
next step, MPB, builds and diagonalizes the Floquet Hamiltonian matrix Hp + Lb and 
calculates the surface amplitudes ( [TTI ) required to construct the i?-matrix (|l^) in the 
inner region. In the current implementation, this stage is parallelized using SCALAPACK 
(Blackford et al 1997). The above steps need only be performed once for fixed laser 
frequency and intensity as they are independent of the collision energy. 

The program MPSCAT solves the scattering problem in the outer region. For each 
collision energy, the i?-matrix from the inner region is transformed into the field- 
dressed channel representation, and from the length into the velocity gauge. Its 
inverse is then propagated from the inner region boundary out to a specified distance 
using a log-derivative propagator (Manolopoulos 1986). The propagated solutions are 
matched with those from the asymptotic expansion ( ^3| ) to give the i^T-matrix, from 
which cross sections can be calculated. 

In the field-free case, the i?-matrix is symmetric by construction. This is also the 
case in the length gauge, but not in the velocity gauge since the Floquet expansion is 
truncated and the transformation matrix A ( p6| ) is not quite unitary. Adopting the 
velocity gauge in the inner region would yield a symmetric i?-matrix but the Buttle 
correction then requires non-diagonal elements that are difficult to evaluate. The 
propagation method we employ however does not suppose the matrix to be symmetric. 
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The transformation from the acceleratfon frame into the velocity gauge is not unitary 
for the same reason. As a consequence, the /^-matrix is not symmetric. In practice, we 
find that this does not afi'ect the cross section, provided the full not quite symmetric 
if-matrix is used to determine the S'-matrix and the contributions of the channels with 
the largest angular momenta or largest (absorption and emission) Floquet indices are 
not included. 

From the computational point of view, calculations for transitions between two 
continuum states arc more demanding than for bound-free transitions since more 
partial waves need to be included. It is also necessary to integrate over large interaction 
distances, especially to obtain converged angular distributions. Many operations can 
be written in matrix form, which can be efficiently performed using the BLAS and 
LAPACK (Anderson et al 1999) libraries. All the calculation parameters, such as the 
number of target states, continuum orbitals, angular momenta, Floquet components 
as well as the propagation distance, can be varied independently to ensure the 
convergence of the results. 

9. Conclusion 

The iZ-matrix Floquet theory provides a powerful tool to study the behaviour of 
a general atomic system in a laser field. The method is fully ab initio and non- 
perturbative and allows in particular an exact treatment of all dynamic and radiative 
couplings. The general approach for studying electron- atom scattering in a laser field 
is similar to the treatment of multiphoton ionization presented in papers I and II, but 
several aspects need more careful consideration, in particular the boundary conditions. 
As we have shown in this paper, it is necessary to define a priori the field-dressed target 
states to identify the collision channels and calculate cross sections. This however has 
the advantage of greatly simplifying the formulae in the outer region. 

As in the standard i?-matrix theory on which our method is based, only one 
electron can be treated in the outer region. The energy of the collisional electron 
and the field strength must not be too high, so that ionization or photoionization 
of the target does not occur. The Floquet approach is furthermore only valid for 
a periodic field, which again limits the field intensity as most powerful lasers are 
pulsed. Extensions of the i?-matrix Floquet theory to treat stronger fields are however 
possible. A time-dependent iZ-matrix theory has already been proposed by Burke and 
Burke (1997) to handle non-periodic fields in a one-dimensional model, which could 
be generalised to three dimensions. There is furthermore no fundamental obstacle in 
combining for instance the i?-matrix pseudo-state method (Bartschat et al 1996) or 
the two-dimensional _R-matrix propagator (Dunseath et al 1996) with the Floquet 
approach to include the double continuum. These methods are computationally 
prohibitive at the moment but could become realistic as the power of computers 
increases. 

The theory presented in this paper has been successfully applied to the study of 
free- free processes and simultaneous electron-photon excitation (Charlo et al 1998, 
Charlo 1999, Terao-Dunseath et al 2001). These were the first accurate non- 
perturbative calculations including a realistic multi-electron representation of the 
collision system. We are currently continuing our study of electron-helium scattering 
in a laser field (CO2, Nd-YAG) and intend to apply the method to other complex 
targets of experimental interest. 
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